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t^^- ■ Abstract 

O 

New precise measurements of D^ and D*^ meson production from B mesons and qq continuum 

^ . events near the T{AS) resonance are presented in this paper. Using the BABAR data recorded in 1999 

D I and 2000 of 20.8 fb~ on-resonance and 2.6 fb~ off-resonance, we measure the inclusive branching 

Q_;_ fractions B{B -^ D+X) = (10.93 ± 0.19 ± 0.58 ± 2.73)% and B{B -^ D*+X) = (7.94 ± 0.82 ± 

1^ ■ 0.72 ± 1.99)%, where the first error is statistical, the second is the systematic error, and the third is 

the error due to the D'^ -^ (p7r~^ branching fraction uncertainty. The branching fractions Z!13{B — > 

d(*)+D(*)) = (5.07 ± 0.09 ± 0.34 ± 1.27)% and UB{B -^ L>*+^W) = (4.07 ± 0.42 ± 0.53 ± 1.02)% 

rS ' have been determined from the measured Di* momentum spectra. 
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1 Introduction 

The study of Dl ' meson production in B decays allows exploration of the mechanisms leading 
to the creation of cS quark pairs. Although several Feynman diagrams could lead to Dg' mesons 
in B decays, the spectator diagram (Fig. |I|) is expected to dominate. In addition, Dl mesons 
could be produced from cc continuum events. It has been pointed out |j|] that the rate from B 
decays may be large. This might help to explain some of the theoretical difficulties ||2| in accounting 
simultaneously for the total inclusive B decay rate and the semileptonic branching fraction of the 
B meson. The measurement of the Dl ' momentum allows a determination of the fraction of 
two-body and multi-body decay modes, which will aid understanding h — > ccs transitions. 

In this paper, measurements of i? — > DfX and B -^ D*^X production rates and momentum 
spectra^ are presented. These mesons are reconstructed using the decays D^ -^ cpir^ and D*+ -^ 




Figure 1: The main spectator diagram leading to the production of Dg mesons in B decays. 



2 The BaBAR detector and data set 

The data used for this analysis were collected with the BABAR detector y] at the PEP-II asymmetric- 
energy collider [ffl] at the Stanford Linear Accelerator center. An integrated luminosity of 20.8 fb~ 
was recorded corresponding to about 22.7 million produced BB pairs at the T(45) resonance ("on- 
resonance") and 2.6 fb~ at an energy about 40MeV below the BB threshold ("off-resonance"). 
Since a detailed description of the BABAR detector is presented in Ref. [^, only the components of 
the detector most crucial to this analysis are briefly summarized below. 

Charged particles are detected and their momenta measured by a combination of a central 
drift chamber (DCH) with a helium-based gas and a five-layer (double-sided) silicon vertex tracker 
(SVT), within a 1.5 T solenoidal field produced by a superconducting magnet. The tracking 
system covers a solid angle of 92% in the center-of-mass frame. Charged particles are identified 
using the ionization energy loss (dB/dx) measured in the DCH and SVT and the Cherenkov 
radiation detected in a ring imaging Cherenkov device (DIRC). Photons are identified by the Csl 
electromagnetic calorimeter. 

^Reference in this paper to a specific decay channel or state also implies the charge conjugated decay or state. 
The notation Ds* means either Dt or Dl^ . B -^ Dg D^*' is a general representation of any of modes with cs, 
cq including their excited states. 



3 Df and Z>*+ selection 

The analysis reported here uses only the decay mode Df — > (J)tt~^ , with (p -^ K^K~ , as this channel 
offers the best signal-to-background ratio. The charged tracks are required to originate from within 
±10 cm along the beam direction and ±1.5 cm in the transverse plane of the interaction point and 
leave at least 12 hits in the drift chamber. 

In order to obtain a sufficiently clean sample, kaon identification is required for the tracks 
forming the (p meson by using dE/dx information from DCH and SVT and the Cherenkov angle 
and the number of photons as measured by the DIRC. The kaon selection is based on the likelihoods 
given by each detector and uses, for each track, the ratio of likelihoods for the pion and the kaon mass 
hypotheses, Lt^/Lk- If this ratio is less than unity for at least one of the considered subsystems, 
the particle is selected as a kaon. The DIRC is used both in a positive identification mode and also 
in a veto mode for the case where a kaon with the measured track momentum would not be above 
the Cherenkov threshold. A tighter level of identification is also available using a total likelihood 
defined as the product of the likelihoods of each subsystem. In this case the track is tagged as a 
kaon if the ratio of the total likelihoods for the pion and kaon mass hypotheses is less than unity. 

Three charged tracks originating from a common vertex are combined to form a D^ candidate. 
Two oppositely charged tracks have to be identified as kaons by satisfying the basic criteria and at 
least one of them has to satisfy the tighter selection. The K^K~ invariant mass must be within 
8MeV/c^ of the nominal cp mass [^. In this particular decay, the (p meson is polarized longitudinally 
and therefore the angular distribution of the kaons has a cos^ 6h dependence, where 6h is the angle 
between the K^ and D^ in the (p rest frame. This angle is required to satisfy |cos0f/| > 0.3, 
thereby keeping 97.5% of the signal while rejecting about 30% of the background. 

Using the selection described above, a clean Df signal of 47794±311 events is observed (Fig. ||). 
A clear signal for the Cabibbo-suppressed decay mode D^ -^ (piT~^ is also observed. 

D*'^ mesons are reconstructed using the decay D*~^ — s- Df^, with the subsequent decay D'^ — > 
(pir^. Df candidates are selected by requiring the (pn invariant mass to be within 2.5 standard 
deviations (a) of the peak value. These are then combined with the "single photons" of the event, 
which are required to satisfy E^ > 50MeV, where E^ is the photon energy in the laboratory frame, 
and E* > 110 MeV, where E* is the photon energy in the T(4S) rest frame. In order to reduce the 
combinatoric background, the candidate photon should not form a vr*^, defined by a total energy 
E* > 200 MeV and an invariant mass 115 < Myy < 155 MeV/ c^, when combined with any other 
photon in the event. The distribution of the mass difference AM = M^+^ — Mr-,+ is shown in 
Fig. |3|. A clear peak with 14392 ± 376 events is observed. 



4 Extraction of the £)(*)+ momentum spectra 

The momentum spectrum of Df mesons in the T(4S') rest frame is extracted by fitting the (pir 
invariant mass distribution in each momentum bin. The bins are chosen to be 200MeV/c wide, 
which is much larger than the momentum resolution (~ 6MeV/c). The fit function is a single 
Gaussian for each of the Df and the D^ signals, with the constraint of a common width. The 
combinatorial background is accounted for by an exponential. As there are many more events for 
the on-resonance data, the number of Df in the off-resonance data is extracted with the same fit 
function but with Mr)+, M^+ and a fixed to the values obtained from the binned chi-squared fit 
to the on-resonance data. 




4000 



M(4)7r), GeV/c 




M(D Y)-M(D ), GeV/c 



Figure 2: The observed (pTr invariant mass 
spectrum. The lower mass peak corresponds 
to the Cabibbo-suppressed decay mode D^ -^ 
(pTT^ . The fit function is a single Gaussian for 
each peak, with their widths constrained to be 
equal, on top of an exponential background. 



Figure 3: Distribution of the AM = M,-.+^ — 
Mp.+ mass difference. The fit function is a 
Crystal Ball function for the signal on top of 
a threshold function, as described in text. 



In the same way as for Df, the momentum spectrum of D*^ mesons in the T{4S) rest frame 
is extracted by fitting the AM invariant mass distribution for 250MeV/c wide momentum bins. 
The AM distribution for the signal is characterized by an asymmetric shape to account for energy 
leakage and shower shape fluctuations. The fit function for the signal is the Crystal Ball function Q. 
For the background, a threshold function 



/(AM) =pi{AM 



■ ■n^V-i ^^{^^i-P2) 



is used, with the four free parameters pi being determined from the fit. After ensuring that the 
connection point between the Gaussian and power-law tail of the Crystal Ball function does not 
depend on momentum and agrees with the Monte Carlo, this parameter has been fixed at 0.89(T 
in the final fit. The off-resonance data are again fit with the signal shape parameters AM and a 
fixed to the values obtained from the fit to the on-resonance data. 

The uncertainty on the shape of the background leads to an additional systematic error. This 
error is estimated by using different parameterizations for the background shape. 

The efficiency obtained from Monte Carlo with generic BB and cc events varies as a function 
of the uf momentum (p*) in the T(4S') rest frame and ranges from 20% when the D^ is at rest 
to 40% for p* = 5GeV/c, and from 5% to 20% for D*'^ . The number of reconstructed D^ and 
D*g^ is corrected bin-by-bin for the efficiency. The efficiency-corrected number of D^ and -D*^ as 
a function of their momentum in the T(4S') rest frame is shown in Fig. 0. 
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Figure 4: The (a) -D^ and (b) D*^ efficiency-corrected momentum spectra for on-resonance data 
(solid circles) and for scaled off-resonance data (open circles). 

5 Branching fractions 

The Df and D*~^ cross sections for production from the qq continuum are obtained by integrating 
the spectrum obtained from the off-resonance data. This gives the preliminary results 



a(e+e~ -^ DfX) x B{D^, 



and 



CT(e+e- -^ D*^X) X B{D 



57r 



OTT 



7.55 ± 0.20 ± 0.34 pb 



5.79 ± 0.66 ± 0.50 pb. 



The off-resonance data are scaled according to the luminosity ratio and then subtracted bin-by- 
bin from the on-resonance data in order to find the Ds momentum spectra from B meson decays. 
It is important to note that, with this method, the result is independent of any assumption about 
the shape of the fragmentation function, and most of the systematic errors due to the background 
parameterization cancel. Integrating the spectrum after continuum subtraction gives a total Df 
yield from B meson decays of 87711ibl485. This corresponds to the inclusive preliminary branching 
fraction 

3.6 ± 0.970 
(10.93 ±0.19 ±0.58) 



BiB ^ DtX) 



B{Dt -^ </'7r+) 



%. 



The total Dl'^ yield from B meson decays is 60047±6201 events, leading to the inclusive preliminary 
branching fraction 



B{B -^ D*+X) 



(7.94 ±0.82 ±0.72) x 



3.6 ±0. 



B{D^s 



t)TT^ 



In the results above, the first error is statistical, the second is the systematic error and the third 
error, which is dominant, is due to the uncertainty in the Df — > (pir'^ branching fraction [^]. The 
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Table 1: Systematic errors for B{B -^ Ds* X) 



Source 


Fractional Error on B (%) 




B^DfX 


B - 


-^D*+X 


Signal shape 


0.5 




3.0 


Background parameterization 


0.4 




4.2 


Monte Carlo statistics 


2.5 




4.2 


Bin width 


1.4 




2.0 


Total for D+ yield 


2.9 




7.0 


Number of BB events 


1.6 




1.6 


B{(I)^K+K~) 


1.6 




1.6 


Particle id efficiency 


1.0 




1.0 


Tracking efficiency 


3.6 




3.6 


BiD*+ ^ D+j) 






2.7 


Photon efficiency 






1.3 


7r° veto 






2.7 


Total systematic error 


5.3 




9.0 



various contributions to the systematic error are listed in Table ||. One of the dominant systematic 
errors is the 3.6% total uncertainty due to our knowledge of the tracking efficiency (1.2% per track 
for the decay chain Df -^ 4>'^~^ , 4> -^ K^K^). 

6 Fits to Z)(*)+ momentuni spectra 

In the T(4S') rest frame, two-body B decays produce Dg mesons with a flat momentum spectrum 
over a 300MeV/c wide range. In B decays, the Dg momentum spectrum is essentially governed 
by the production of direct Dg . Other cs states such as -Dsi(2536) and D*2(2573) primarily 
decay to D^*'K. Because D*~^ decays to Dfj or Df-K^, the Df momentum distribution is slightly 

broader and shifted downward compared to direct production from B — > D^X. 

(*)+ 
In fitting the observed momentum spectra, three different sources of Ds mesons in B decays 

are considered: 

(1) B — > Ds D^*> decays. The relative branching fractions of the individual channels can be 
taken either from existing measurements Q or from predictions assuming factorization ||8|, |9|, 
IC]. The fit is performed for both cases, with the assumption /£)•+ = /£>+ for the theoretical 



models, where /„(.)+ are the Dl decay constants. 

(2) B -^ Di*'>^D** decays. The contributions from B -^ L»i*^+D^(j = 1/2), B -^ L>i*^'^^i(2420), 
B -^ Di*^^Di{j = 1/2) and B -^ 1)1*^+^^(2460) are included in this source. 

(3) Three-body B -^ Dg* D^*'tt/p/u! decays. Since little is known on these decays, they are 
attributed an equal weight and the momentum distributions are generated according to phase 
space. 
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Figure 5: The fit result for (a) Df and (b) Ds'^ momentum spectra. The data are dots with 
error bars, the histograms are the components of the fit function described in the text. Type (1) 
is 5 ^ Z)i*)+^(*), Type (2) is S ^ D^*^^D** and Type (3) is 5 ^ D^*^^ D^*\ / p / oj . The soHd 
histogram is the sum of three components. 



As a result of a chi-squared fit of Df momentum spectrum with these components, the ratio of 
two-body modes to the total inclusive rate is determined to be 



^ -^—^ — '- = 46.4 ± 1.9 ± 0.6 %, 



where the first error is statistical and the second is due to model uncertainty. This last is obtained 
from the variation of the fit result with different individual contributions from the modes included 
with each of the three sources of Dj" mesons, as discussed below. From the fit to the D\^ momentum 
spectrum, we find 

where the first error is statistical, the second error represents the systematic uncertainty due to the 
background parameterization (negligible for B —>■ D), D^*>), and the third error is due to model 
uncertainty obtained as for B -^ Dg* D^*' . 

The fit is performed under different assumptions for the relative contributions of the modes 
in source (1), varied according to the theoretical predictions and measurements. Different weights 
of i? ^ D^D** and B -^ D*^D**, as well as different relative branching fractions of the four 
modes within source (2), are also used. In source (3), two cases are considered: either B — > 
Ds* D^*>7r oi B ^ Ds D^*'p/uj is assumed to be dominant. The best x^ foi' the fit to the inclusive 
D*^ momentum spectrum is obtained when the contribution from B -^ Dg* D^*' p/to is dominant 
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compared to B ^ Dg D^*'7r. The results of the fits to the Ds momentum spectra are shown in 
Fig. |5| for one of the assumptions. 

Using the fit results and the relative rates for B -^ Ds* D^*' we find the preliminary results 

UB{B -^ Z)W+^{*)) = (5.07 ± 0.09 ± 0.34 ± 1.27)%, 

UB{B -^ Dl+D'^*^) = (4.07 ± 0.42 ± 0.53 ± 1.02)%. 

where the errors from the fits to the momentum spectra are added in quadrature with the systematic 
error due to the Df -^ (pir'^ branching fraction uncertainty. 

7 Summary 

In summary, preliminary branching fractions for inclusive B ^ Ds X production have been 
determined to be B{B -^ D+X) = (10.93 ± 0.19 ± 0.58 ± 2.73)% and B[B -^ D*+X) = (7.94 ± 
0.82 lb 0.72 zb 1.99)%. The Ds cross sections from qq continuum events at about 40MeV below 
r(4S') mass are ^(e+e- -^ DfX) x B{D+ -^ 4>tt+) = 7.55 ± 0.20 ± 0.34 pb and a{e+e- -^ 
D*'^X) X B{Df -^ (/)7r+) = 5.79 ± 0.66 ± 0.50 pb. Our results for Z?+ are in agreement with 
previous measurements R, |ll|, although with considerable improvement in accuracy. In contrast 
to previous measurements, our results do not rely on any assumptions concerning the shape of the 
fragmentation function. 

From a fit to the Dl ' momentum spectra, preliminary results have been obtained for the 
fraction of all two-body B -^ Ds* D^*' decays relative to the total inclusive Df yield (46.4 it 1.9 ± 
0.6)% and for all B -^ D*+D^*'^ decays relative to the total inclusive D*+ yield (53.3 ± 4.5 ± 1.6 ± 
2.1)%, where the last error includes the model dependence. Combining these results gives EB{B — > 
^(*)+^(*)) ^ (5 Q7 ^ 0.09 ± 0.34 ± 1.27)% and EB{B -^ Z)*+dW) = (4.07 ± 0.42 ± 0.53 ± 1.02)% 
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